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ABSTRACT 



Recent calculations are presented of top quark polarization in tt pair production 
close to threshold. S-P-wave interference gives contributions to all components of 
fSI \ the top quark polarization vector. Rescattering of the decay products is considered. 

■ Moments of the fourmomentum of the charged lepton in semileptonic top decays 

are calculated and shown to be very sensitive to the top quark polarization. 

o 



1. Introduction 

Threshold production of top quarks at a future electron-positron collider will al- 
low to study their properties with extremely high precision. The dynamics of the 
top quark is strongly influenced by its large width r t w 1.5 GeV. Individual quarko- 



nium resonances can no longer be resolved, hadronization effects are irrelevant and 
an effective cutoff of the large distance (small momentum) part of the hadronic in- 
teraction is introduced™™. This in turn allows to measure the short distance part 
of the potential, leading to a precise determination of the strong coupling constants. 
The analysis of the total cross section combined with its momentum distribution will 
determine its mass with an accuracy of at least 300 MeV and its width to about 10%. 
For a Higgs boson mass of order 100 GeV even the ttH Yukawa coupling could be 
indirectly deduced from its contribution to the vertex correction^. Additional con- 
straints on these parameters can be derived from the forward-backward asymmetry 
of top quarks and from measurements of the top quark spin. Close to threshold, for 
E = y/s — 2m t <C m t , the total cross section and similarly the momentum distribu- 
tion of the quarks are essentially governed by the S'-wave amplitude, with P-waves 



suppressed ~ (3 2 ~ \JE 2 + T 2 /m t ~ 10~ 2 . The forward-backward asymmetry and, 
likewise, the transverse component of the top quark spin originate from the interfer- 
ence between S- and P-wave amplitudes and are, therefore, of order j3 ~ 10 _1 even 
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close to threshold. Note that the expectation value of the momentum is always differ- 
ent from zero as a consequence of the large top width and the uncertainty principle, 
even for E = 0. 

It has been demonstrated™ that the Green function technique is particularly suited 
to calculate the total cross section in the threshold region. The method has been 
extended™ to predict the top quark momentum distribution. A further generalization 
then leads to the inclusion of P-waves and, as a consequence, allows to predict the 
forward-backward asymmetry! It has been shownHthat the same function ip R (p,E) 
which results from the 5-P-wave interference governs the dynamical behaviour of the 
forward-backward asymmetry as well as the angular dependence of the transverse 
part of the top quark polarization. The close relation between this result and the 
tree level prediction, expanded up to linear terms in (3, has been emphasised. The 
relative importances of Z versus 7 and of axial versus vector couplings depend on 
the electron (and/or positron) beam polarization. All predictions can, therefore, 
be further tested by exploiting their dependence on beam polarization. In fact the 
reaction e + e~ — » ti with longitudinally polarized beams is the most efficient and 
flexible source of polarized top quarks. At the same time the longitudinal polarization 
of the electron beam is an obvious option for a future linear collider. Recently^ these 
results have been expanded in two directions: 

• Normal polarization. Calculation of the polarization normal to the production 
plane is a straightforward extension of the previous workfa and is based on the 
same nonrelativistic Green function as before, involving, however, the imaginary 
part of the interference term ifj(p, E). A component of the top quark polariza- 
tion normal to the production plane may also be induced by time reversal odd 
components of the ytt- or Ztt-coupling with an electric dipole moment as most 
prominent example. Such an effect would be a clear signal for physics beyond 
the standard model. The relative sign of particle versus antiparticle polariza- 
tions is opposite for the QCD-induced and the T-odd terms respectively, which 
allows to discriminate the two effects. Nevertheless it is clear that a complete 
understanding of the QCD-induced component is mandatory for a convincing 
analysis of the T-odd contribution. 

• Rescattering. Both t quark and t antiquark are unstable and decay into W + b 
and W~b, respectively. Neither b nor b can be considered as freely propagating 
particles. Rescattering in the tb and bi systems affects not only the momenta 
of the decay products but also the polarization of the top quark. Moreover, in 
the latter case, when the top quark decays first and its colored decay product 
b is rescattered in a Coulomb-like chromostatic potential of the spectator t, 
the top polarization is not a well defined quantity. Instead one can consider 
other quantities, like the total angular momentum of the Wb subsystem, which 
are equal to the spin of top quark in the situation when rescattering is absent. 
These rescattering corrections are suppressed by a s . The resulting modifications 
of the momentum distribution are therefore relatively minor and as far as the 
total cross section is concerned can even be shown to vanis h§. In contrast the 
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forward-backward asymmetry as well as the transverse and normal parts of the 
top quark spin are suppressed by a factor ~ j3. Thus, they are relatively more 
sensitive towards rescattering corrections. Rescattering in the bb system is less 
important and will be neglected. 

It is well known0 that the direction of the charged lepton in semiieptonic decays 
is the best polarization analyzer for the top quark. The reason ig!2l that in the top 
quark rest frame the double differential energy-angular distribution of the charged 
lepton is a product of the energy and the angular dependent factors. The angular 
dependence is of the form (1+P cos 6), where P denotes the top quark polarization and 
9 is the angle between the polarization three-vector and the direction of the charged 
lepton. Gluon radiation and virtual corrections in the top quark decay practically 
do not affect these welcome properties^. It is therefore quite natural to perform 
polarization studies by measuring the inclusive distributions of say fi + in the process 
e + e~ — > t{^ + v^b)i{jets). This can be also convenient from the experimental point 
of view because there is no missing energy-momentum for the i subsystem. From 
the theoretical point of view the direction of the charged lepton can be considered as 
another quantity which is equivalent to the top quark polarization when rescattering 
is absent. Of course, it is well defined also in the case of bt rescattering. However, 
the semi-analytic calculation of the latter contribution is a .very difficult task because 
production and decay mechanisms are coupled. A way outt9 is to calculate moments 
of the charged lepton four-momentum distributions. The results of this analysis are 
published elsewhere^. 

2. Green functions, angular distributions and quark polarization 

2.1. The nonrelativistic limit 

Top quark production in the threshold region is conveniently described by the 
Green function method which allows to introduce in a natural way the effects of the 
large top decay rate T t and avoids the summation of many overlapping resonances. 
The total cross section can be obtained from the imaginary part of the Green function 
G(x = 0,x' = 0, E) via the optical theorem. To predict the differential momentum 
distribution, however, the complete x-dependence of G(x, x' = 0,E) (or, more pre- 
cisely, its Fourier transformed) is required. In a calculation with non-interacting 
quarks close to threshold the forward-backward asymmetry, the leading angular de- 
pendent term ~ cos$ and the transverse part of the top quark polarization are all 
proportional to the quark velocity f3 and originate from the interference of a S'-wave 
with a P-wave amplitude. These distributions are described by V'-G(x, x', E)\ x > =0 or, 
equivalently, by the component of the Green function with angular momentum one. 
The connection between the relativistic treatment and the nonrelativistic Lippmann- 
Schwinger equation has been discussed in the literaturepfl. The subsequent discussion 
follows these lines. It includes, however, also the spin degrees of freedom and is, fur- 
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Figure 1: Lippmann-Schwinger equation in diagrammatical form. 



thermore, formulated sufficiently general such that it is immediately applicable to 
other reactions of interest. The main ingredient in the derivation of the nonrelativis- 
tic limit is the ladder approximation for the vertex function r^. This vertex function 
is the solution of the following integral equation depicted in Fig.Q: 

T c = C + J-0^ (-I 4 ™*) D ^ p - *)7"£f(* + §)Ic(*, q)S F (k - | )-f, (1) 
with C = or 7^75 in the cases of interest for e + e~-annihilation. The conventions 



q/2 + p 



k + q/2. 
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Figure 2: Definition of the four- momenta. 



for the flow of momenta are illustrated in Fig|| The four-momenta are related to the 
nonrelativistic variables by 

q = (2m t + £,0) 

P = (0,p) (2) 
k = (k ,k). 

In perturbation theory the ladder approximation is motivated by the observation that 
for each additional rung the energy denominator after loop integration compensates 
the coupling constant attached to the gluon propagator. This is demonstrated most 
easily in Coulomb gauge. Contributions from transverse gluons as well as those from 
other diagrams are suppressed by higher powers of (3 ~ a s . The gluon propagator is 
thus replaced by the instantaneous nonrelativistic potential 

4 

- A-nasD^p) -> i V(p) 5 M o<W (3) 
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The dominant contribution to the integral originates from the region where |k| <C m t . 
Including terms linear in k, quark and antiquark propagators are approximated by 



S F (k + 



T " E 4- h k2 -1- i Tt 



2' #-A-o-^: + / 

A 



2 u 2m t 1 2 

1± 7 ° 



The "elementary" vertex C is independent of ho- (Within the present approximations 
this is even true if C does depend on k as is the case in the analogous treatment of 
77 — > ti discussed below.) Up to and including order (3 terms a selfconsistent solution 
of the integral equation (ffF) can be obtained if Yq is taken independent of k and the 
nonrelativistic spins of t and i. The ko integration is then easily performed and the 
integral equation simplified to 



^V(p - k) 7 °(A + - ^) Tc ^ E) (A, - ^) 7 °. 



nit 



In the calculation of the cross section for the production of t plus t with momenta 
q/2 ±p and spins s± respectively traces TC of the following structure will arise: 

H = Tr{V + (^+ P ,s + )TcV-(^-p,s^tc}, (6) 

V± ^ S) = 1^ 2—' (7) 

where we allowed for mixed terms with C different from C . Expanding again up to 
terms linear in k, this trace can be transformed into 

H = Tr{<S + f c <S_iV}, (8) 
1 ± s± ■ £ 

S± = | , (9) 

£ = 7757°= (o °), (10) 
with the nonrelativistic reduction defined through 

r c (p^) = A + (i-B_2)r c( p^)(i-H_l) A _. 

It is thus sufficient to calculate the "reduced" vertex function IV Dropping again 
terms of order k 2 , the corresponding integral equation is cast into a particularly simple 
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form 



r c (p,£)=c( P ) + 



d 3 k f c (k,E) 



E-^ + iT t 



(12) 



Consistent with the nonrelativistic approximation only the constant and the linear 
term in the Taylor expansion of the elementary vertex will be considered* 



C(p)=C(0) + D-p 



(13) 



The matrices C(0) and D may in general depend on external momenta, polarization 
vectors or Lorentz indices. A selfconsistent solution for the vertex Tc is then given 
by 

r c (p, E) = C(0)/C s (p, E) + D ■ p £ P (p, E ) (14) 
The scalar vertex functions /Cs,p depend on 

p = IpI 

and E only. They are solutions of the nonrelativistic integral equations 
Ks(p,E) = 1 + 



d 3 k )C s (k,E) 
d 3 k p • k 



(2tt) 3 p 2 



E-^ + ir t 

mt 1 

V(p - k) ^ 



(15) 
(16) 



and are closely related to the Green function Q(p, x, E) which, in turn, is a solution 
of the Lippmann-Schwinger equation 



E-^ + iT t 

m t 



£(p,x,£) = e* px + 



d 3 k 



V(p-k)£(k,x,£). 



;i7) 



Let us denote the first two terms of the Taylor series with respect to x by G and F 
respectively: 

0(p, x, E) = G(p, E) + x ■ p F(p, E) + ... 
They are solutions of the integral equations 



f d h 

G(p,E) = G (p,E) + G (p,E)J—-V(p-k)G(k,E) 



with 



F(p,E) = G (p,E) + G (p,E) 



G (p,E) 



d 3 k 

d 3 k p • k 

(2tt) 3 p 2 

1 



V(p-k)F(k,E), 



E-^ + iTt 



In the notation of Kiihn et al.0 one gets: C(0) = A + o A_ and D = A+ [--L-{0 O , 7}+ + O 



(18) 
(19) 

(20) 
A_. 
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and the relation between Green function and vertex function 

G(p, E) = G (-p, E)JC s (j>, E), F(p, E) = G (p, £)Mp, E) (21) 

is evident. In the case of e + e~-annihilation top production proceeds through the 
space components of the vector and axial vector current. The relevant elementary 
vertex C(p) is given by 

7KP) = A +7 ,-A_ (22) 
77y 5 ( P ) = M^)(7xp)jA- (23) 

for vector and axial current respectively. Production of ti in 77-fusion would lead to 
an elementary vertex of the form 

C(0) oc z(ei x e 2 )n e -A +75 A_ (24) 

D oc — A+ [(ei ■ e 2 )(n e - • j)n e - + (e 2 • 7) e 1 + (e x ■ 7) e 2 ] A-, (25) 
m t 

with 61, e 2 the polarization vectors of the photons, and the present formalism applies 
equally well. This case has been studied by Fadin et al.li3. 

2.2. Top production in electron positron annihilation 

With these ingredients it is straightforward to calculate the differential momen- 
tum distribution and the polarization of top quarks produced in electron positron 
annihilation. Let us introduce the following conventions for the fermion couplings 

v f = 21) - Aq f sin 2 6 W , a f = 2I 3 f . (26) 

P± denotes the longitudinal electron/positron polarization and 

*=£i£ (27) 

can be interpreted as effective longitudinal polarization of the virtual intermediate 
photon or Z boson. The following abbreviations will be useful below: 

oi = QeQt + ( v e + a l) v td 2 + 2q e q t v e v t d 
a 2 = 2v e a e v 2 t d 2 + 2q e q t a e v t d 
a 3 = 4v e a e v t a t d 2 + 2q e q t a e a t d 
a 4 = 2(f 2 + a 2 e )v t a t d 2 + 2q e q t v e a t d 
1 s 



d 



16 sin 2 9 W cos 2 9 W s - Ml 
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The differential cross section, summed over polarizations of quarks and expanded up 
to terms linear in p, is thus given by 



dp Awm 



2r / 16a 

a 1 + X a 2 )(l-—^)\G(p,E)\< + 



da 3a2r S-p + p- 



3tt 



+ (a 3 + x« 4 ) ( 1 - ^) —Re ( G(p, E)F*(p, E) ) cos 
Sir J m t 



(28) 

The vertex corrections from hard gluon exchange for 5*-waveli3 and P-wave@ am- 
plitudes are included in this formula. It leads to the following forward-backward 
asymmetry 

A FB (P,E) = C FB ( X )V R (V,E), (29) 

with 

1 a 3 + x«4 



tp R (p,E) = Reip, and 



Cfb(x) = -^-^, (30) 
2 ai + X°2 



_ (1-W37T) p F*(p, g 
(1 - 8a s /37r) m t G*(p, 



This result is still differential in the top quark momentum. Replacing (p(p, E) by 

= (1 - 4a,/37r) /o P ^pg^(P^Mp^) 
1 J (l-8a,/37r) J P ™ c^pp 2 |G(p, E)\ 2 ^ 1 J 

one obtains the integrated forward-backward asymmetry! The cutoff p m must be 
introduced to eliminate the logarithmic divergence of the integral. For free particles 
(or sufficiently far above threshold) one finds for example 

= J^- + ^2E^m Pm /m t (33) 
V m t \Jm t Eix 

This logarithmic divergence is a consequence of the fact that the nonrelativistic ap- 
proximation is used outside its range of validity. On may either choose a cutoff of order 



m t or replace the nonrelativistic phase space element ~pdp/m t by pdp/ym 2 + p 2 . In 
practical applications a cutoff will be introduced by the experimental procedure used 
to define tt-e vents. 



^For the case without beam polarization this coincides with the earlier resultB, as far as the Green 
function is concerned. It differs, however, in the correction originating from hard gluon exchange. 



s 



2.3. Polarization 



To describe top quark production in the threshold region it is convenient to align 
the reference system with the beam direction (Fig|J) and to define 



n e - 



x n, 



(34) 



sn x sy. 

In the limit of small j3 the quark spin is essentially aligned with the beam direction 




Figure 3: Definition of the spin 
directions. The normal compo- 
nent sn points out of the plane. 



apart from small corrections proportional to (3, which depend on the production angle. 
A system of reference with sy defined with respect to the top quark momentum!^ is 
convenient in the high energy limit but evidently becomes less convenient close to 
threshold. 

Including the QCD potential one obtains for the three components of the polarization 



V\\{v,E, X ) = Cf( X ) + Cf(x)^H(P>^) costf 
V ± (p,E, X ) = C ± (x)Vn(P,E)sin$ 
P N (jp,E,x) = C N (xV 7 (p,£)sintf, 



(35) 
(36) 
(37) 



cj(x) = -*L±^, C\{ X ) = (i - x 2 ) ; 2Q3 ~ Ql ^ , 

ri t \ Ifl4 + X«3 n f x la 3 + xa 4 , . , ,N | 

C±(x) = -~ — — , C N (x) = -- — — = -Cfb(x), 

2 ai + x«2 2 ai + xa 2 

with <^(p, -E) = Im<^, and ip(p,E) is defined in (|31|). The momentum integrated 
quantities are obtained by the replacement <p(p,E) — ► $(£"). The case of non- 
interacting stable quarks is recovered by the replacement $ — > /3, an obvious con- 
sequence of fl3~2"|). 

Let us emphasize the main qualitative features of the result. 

• Top quarks in the threshold region are highly polarized. Even for unpolarized 
beams the longitudinal polarization amounts to about —0.41 and reaches ±1 for 
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fully polarized electron beams. This later feature is of purely kinematical origin 
and independent of the structure of top quark couplings. Precision studies of 
polarized top decays are therefore feasible. 

• Corrections to this idealized picture arise from the small admixture of P-waves. 
The transverse and the normal components of the polarization are of order 10%. 
The angular dependent part of the parallel polarization is even more suppressed. 
Moreover, as a consequence of the angular dependence its contribution vanishes 
upon angular integration. 

• The QCD dynamics is solely contained in the functions ip or $ which is the same 
for the angular distribution and the various components of the polarization. 
However, this "universality" is affected by the rescattering correctionsl!2l. These 
functions which evidently depend on QCD dynamics can thus be studied in a 
variety of ways. 

• The relative importance of P-waves increases with energy, $ ~ yj ' E /m t . This 
is expected from the close analogy between Re$ and (3. 

The Ci are displayed in Fig.[| as functions of the polarization x- For the weak mixing 
angle a value sm 2 9w = 0.2317 is adopted, for the top mass m t = 180 GeV. As 
discussed before, CP assumes its maximal value ±1 for x — Tl an d the coefficient Cj 
is small throughout. The coefficient C± varies between +0.7 and —0.5 whereas Cn is 
typically around —0.5. The dynamical factors $ are around 0.1 or larger, such that 
the P-wave induced effects should be observable experimentally. 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 



X X 

Figure 4: The coefficients © for = m t = 180 GeV. 

The normal component of the polarization which is proportional to tpj has been pre- 
dicted for stable quarks in the framework of perturbative QCD HE. In the threshold 
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region the phase can be traced to the ti rescattering by the QCD potential. For a 
pure Coulomb potential V = —4a s /3r and stable quarks the nonrelativistic problem 
can be solved analytically li^l and one finds 

The component of the polarization normal to the production plane is thus approxi- 
mately independent of E and essentially measures the strong coupling constant. In 
fact one can argue that this is a unique way to get a handle on the scattering of heavy 
quarks through the QCD potential. 
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